Distribution of local density of states in superstatistical random 

matrix theory 

A. Y. Abul-Magd 

Department of Mathematics, Faculty of Science, Zagazig University, Zagazig, Egypt 

(Date textFebruary 6, 2008) 

Abstract 

We expose an interesting connection between the distribution of local spectral density of states 
arising in the theory of disordered systems and the notion of superstatistics introduced by Beck 
and Cohen and recently incorporated in random matrix theory. The latter represents the matrix- 
element joint probability density function as an average of the corresponding quantity in the stan- 
dard random-matrix theory over a distribution of level densities. We show that this distribution 
is in reasonable agreement with the numerical calculation for a disordered wire, which suggests to 
use the results of theory of disordered conductors in estimating the parameter distribution of the 
superstatistical random-matrix ensemble. 
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I. INTRODUCTION 



The formalism of superstatistics (statistics of a statistics), has recently been proposed by 
Beck and Cohen |l|] as a possible generalization of statistical mechanics. Superstatistics arises 
as weighted averages of ordinary statistics (the Boltzmann factor) due to fluctuations of one 
or more intensive parameter (e.g. the inverse temperature). It considers a non-equilibrium 
system as traveling within its phase space which is partitioned into cells. Within each 
cell, the system is described by ordinary Maxwell-Boltzmann statistical mechanics, i.e., its 
statistical distribution is the canonical one e~ f3E , but (3 varies from cell to cell, with its own 
probability density f(f3). This formalism has been elaborated and applied successfully to a 
wide variety of physical problems, e.g., in 0, [], [5, [fi, 7], S El Q, 



Superstatistics has been applied to model systems with partially chaotic classical dynam- 
ics within the framework of random-matrix theory (RMT) in Ref. [12I. \\\. It has provides 
a possible mechanism for the initial stage of transition of a system out of the state of chaos 
but fails to reproduce the Poisson statistics that is believed to describe regular systems. In 



the standard RMT 
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171 ]. a chaotic system is modeled by an ensemble of random 



matrices that depends only on the symmetry of the system. For example, a system of spin- 
less particles, which has a time-reversal symmetry is represented by a Gaussian orthogonal 
ensemble (GOE). The joint distribution of matrix elements of the Hamiltonian H is pro- 
portional to exp [— 77 Tr (if^iJ)] , where Tr is trace and W is the Hermitian conjugate of H. 
The parameter rj is related to the square of the mean level density. This distribution is based 
on two main assumptions: (i) the matrix elements are independent identically-distributed 
random variables, and (ii) their distribution is invariant under unitary transformations. For 
most of the physical systems, however, the phase space is partitioned into regular and chaotic 
domains. These systems are known as mixed systems. Attempts to generalize RMT to de- 



16| . Most of these attempts are 



scribe such mixed systems are numerous; for a review see 
based on constructing ensembles of random matrices whose elements are independent but 
not identically distributed, e.g. in Q, Q, 20, 21]. Thus, the resulting expressions are 
not invariant under base transformation. The superstatistical generalization follows another 
route. It keeps base invariance, but violates matrix-elements independence. The intuitive 
explanation for using superstatistics is based on the ansatz that the spectrum of the sys- 
tem under consideration is partitioned into small cells. Within each cell, the spectrum 
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is described by an ordinary Gaussian random-matrix ensemble, but ensemble parameter 77 
varies from cell to cell.. One may define the density of states constituting a single cell as 
a 'local density of states' (LDOS) proportional toy 7 //. Superstatistics assumes that LDOS 
and thus the parameter 77 is no more a constant parameter as in the original RMT, but 
allowed to fluctuate according to a distribution f(rj). The joint matrix-element distribution 
is represented as an average over exp [— 77 Tr [H^HX\ with respect to the parameter 77. The 
resulting distribution depends on the matrix elements in the form of Tr( H^Hj, which is 
base independent as the corresponding distribution in ordinary RMT. 

The central question for superstatistics is how to choose the parameter distribution. 
In superstatistical thermodynamics, one obtains Tsallis' statistics 2J, |25[ when the inverse 
temperature (3 has a x 2 distribution, but this is not the only possible choice. Beck and Cohen 
give several other possible examples of functions which are possible candidates for /(/?). 
Generalized entropies, which are analogous to the Tsallis entropy, can be defined for these 
general superstatistics 261. l27l|. Sattin 6] suggested that, lacking any further information, 



the most probable realization of / 



' 0) will be the one that maximizes the Shannon entropy. 



This latter approach was used in [13j to estimate the distribution f{rf) of the parameter 77 
of RMT. 

The object of the present paper is to find out whether the local density of states defined 
above is related to the local spectral density of states (in the literature, also LDOS) which is 
relevant to many practical applications in the field of condensed matter physics 28] . There, 
the notion of LDOS follows from the conjecture that even in a metallic sample there is a finite 
probability to find "almost localized" eigenstates so that the density of eigenstates is different 
at different locations in the sample. Here, by localization we shall mean a situation, in which 
eigenfunctions are localized in the space of "unperturbed" eigenvalues on a scale which is 
significantly smaller than the size of the ensemble. Indeed, the size of the region which is 
populated by an eigenfunction (termed localization length) measures the maximum number 
of basis state coupled by perturbation (off-diagonal elements). In the state of chaos, the 
localization length approaches the size of the ensemble, which means that the eigenfunctions 
become ergodic, i.e., extended over the whole energy shell. Then, the LDOS coincides with 
the "global" density of states and both follow Wigner's semicircular law. As the system 
departs from chaos, the localization length decreases and the LDOS becomes distinguished 
from the global one. The main question, which we seek to answer is the following: "Can 
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we benefit from the achievements of the well developed theory of disordered conductors 
in estimating the parameter distribution f(rj) of the superstatistical ensemble?" Instead 
of answering this question directly (which is technically quite difficult), we compare the 
superstatistical distribution f(rj) derived in [l]] from the principle of maximum entropy with 
the distribution of LDOS deduced by Altshuler and Prigodin [29] for disordered conductors. 
After a brief review of our previous results, given in Section II, we show in Section III that 
the two distributions have similar shapes especially near the chaotic limit. The distribution 
f(r]) is com par ed with the distribution of LDOS obtained in the numerical simulation of the 
closed wire |30| . While the agreement of the superstatistical distribution with the numerical 
experiment is not completely satisfactory, it still demonstrates the analogy between the 
notions of LDOS in the two disciplines. Section IV shows by comparison with a numerical 
experiment [31] that the level-density distribution obtained for a random wire can be used 
as a parameter distribution in the superstatistical random matrix theory. The conclusion of 
this work is formulated in Section V. 



II. FORMALISM 



For the sake of completeness and clarity, we recall the derivation of the superstatistical 



model introduced in [13]. RMT models the Hamiltonian of a chaotic system in terms of an 
ensemble of random matrices, whose matrix elements have a Gaussian probability density 
distribution 

P (H) oc exp [-r/Tr (H^H)] . (1) 
The mean level density for a GOE with a large dimension N is given by Wigner's semi-circle 



law [14 1 



p( v ,e) = ±-[Nr ] (l-r ] e 2 /N)} 1/2 6(1 - V e 2 /N), (2) 

where e is the eigenvalue of H and 0(A) is the Heaviside step function. In practical calcu- 
lations with GOE, one usually avoids the ends of the spectrum and works in a region with 
a nearly constant level density equal to p(rj, 0) ~ ^/rj. 

The superstatistical generalization models the quantum-number space of a system of 
mixed regular- chaotic dynamics as made up of many smaller cells that are temporarily in 
a chaotic phase. Each cell is large enough to obey the statistical requirements of RMT but 
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has a different distribution parameter r\ associated with it, according to a probability den- 
sity f(rj). Consequently, the superstatistical random-matrix ensemble describes the mixed 
system as a mixture of Gaussian ensembles. Its matrix-element joint probability density 
distributions obtained by integrating distributions of the form in Eq. (1) over all positive 
values of 77 with a statistical weight f{r)), 

pw=f /(, r [ -y g)] *,. <3) 

where Z{rj) = J exp [— i]Tr (H^H)] drj. Here we use the "B-type superstatistics" [lj. The 
parameter 77 may be expressed in terms of the local mean level spacing D as 

D = ^ (4) 

where c is a constant depending on the size of the ensemble, which can be evaluated by 
setting D = l/p(rj,0) and using Eq. (2). 

In the new framework of RMT provided by superstatistics, the local mean spacing D is 
no longer a fixed parameter but it is a stochastic variable with probability distribution f(D). 
Instead, the observed mean level spacing is just its expectation value. The fluctuation of 
the local mean spacing is due to the correlation of the matrix elements which disappears 
for chaotic systems. In the absence of these fluctuations, f(D) = 5{D — 1) and we obtain 
the standard RMT. Within the superstatistics framework, we can express any statistic a 
of a (sufficiently chaotic) mixed system that can in principle be obtained from the joint 
eigenvalue distribution by integration over some of the eigenvalues, in terms of the corre- 
sponding statistic cr^ G '(D) for a Gaussian random ensemble with mean level spacing D. The 
superstatistical generalization is given by 



a 



f(D)a^ G \D)dD. (5) 



The remaining task of superstatistics is the computation of the distribution f(D). Following 
Sattin we use the principle of maximum entropy (MaxEnt) to evaluate the distribution 
f{D). Lacking a detailed information about the mechanism causing the deviation from the 
prediction of RMT, the most probable realization of f(D) will be the one that extremizes 
the Shannon entropy 

POO 

S = - f{D)\nf{D)dD (6) 
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with the following constraints: (i) The major parameter of RMT is r\ defined in Eq. (1). Su- 
perstatistics was introduced in Eq. (3) by allowing 77 to fluctuate around a fixed mean value 
(77). This requires the existence of the mean inverse square of D, (D~ 2 ) = j °° f(D)D~ 2 dD. 
(ii) The fluctuation properties are usually defined for unfolded spectra, which have a unit 
mean level spacing. We thus require J °° f(D)DdD = 1. As a result, we obtain 



/MaxEnt(-D) = C exp 



[2D 

■ tt U 



D 2 
D 2 



(7) 



where a and Dq are parameters, which can be expressed in terms of the Lagrange multipliers 
of the constrained extremization, and C is a normalization constant. Substituting this 
distribution into Eq. (5) one obtains expressions for the (global) level density, the nearest- 
neighbor-spacing (NNS) distribution and the two- level correlation function 13j. Eq. (5) 
has recently been applied 34j to obtain a generalization of the well-known Porter-Thomas 
distribution of transition intensities, relevant for chaotic regimes, for systems with mixed 
regular-chaotic dynamics. This generalization agrees with the data better than currently 
available results that can not explain observed shift of the peak position as the system 
evolves out of the state of chaos. 



III. DISTRIBUTION OF LOCAL DENSITY OF STATES 



The local density of states plays a central role in the theory of disordered metals 28j, |35 |. 
It also known in nuclear physics as strength function [361] . It gives the distribution of basis 
eigenfunctions in terms of eigenstates of the system. It is obtained by projecting a basis 
state k onto exact eigenstates i of the Hamiltonian H and then defined in terms of the 
expansion coefficients C\ as 

v k (E)=]cjfp(E), (8) 

with p(E) as the density of exact eigenstates. Here, the average is taken over a small 
window of the eigenstates i with energies around E. The LDOS has a well-defined classical 



interpretation as shown by Benet et al. [391 ] . The unperturbed energy Eq is not constant 
along a classical trajectory of the full Hamiltonian with a given total energy H = E. If we 
keep the unperturbed energy E fixed, the bundle of trajectories of the total Hamiltonian 
H, which reach the surface of the unperturbed Hamiltonian Hq = Eq, has a distribution 
in the total energy E which is described by a measure ue (E). In the quantum case, this 



measure corresponds to the imaginary part of the retarded Green's function at energy Eq 40]. 
Distributions of LDOS are relevant for description of fluctuations of tunneling conductance 
across a quantum dot |37| . of transport phenomena in disordered wires [38(, of properties 
of some atomic spectra [41(, and of a shape of NMR line 4^]. In random media the LDOS 
is a random quantity. Altshuler and Prigodin [29] have studied the LDOS distribution 
in strictly one-dimensional disordered chains. They obtain closed-form expressions for the 
LDOS distributions in open and closed wires. The normalized distribution for the closed 
system reads 



00 V 
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(9) 
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where oo is a positive parameter. The distribution Pap{uj, uap] v) has a unit mean value 
z/ap and a variance equal to l/oo. We note that this distribution fulfils the two conditions 
imposed by Beck and Cohen l| on the parameter distribution in superstatistics, which are 
normalization and tending to a delta function as oo — > oo. 

Let us now assume that the LDOS v is equal to the inverse of local mean level spacing D. 
In order to justify this assumption, we compare the distribution Pap{v) with the distribution 



MaxEn 



RMT 



[DYo 
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the local mean level densities (7), which has been applied in the superstatistical 
34|. For this purpose, we have to translate the distribution /MaxEnt(-D) obtained 
in the previous section into the context of the problem of LDOS in complex systems. This 
can be straightforwardly done by substituting v — 1/D in Eq. (7). The distribution of 



LDOS in the superstatistical RMT is then given by 



1 



PviaxEntO) = C — exp 



1^ 



-a: 



2z/ 



v 



/ J 



where u = 1/D Q . The normalization constant is given by 



C = 



a' 



(10) 



(11) 



where Gq^ (z \b\, b 2 , 63) is Meijer's G-function [32 
that the average value of the LDOS is equal to 1, which yields 

Gi(« 3 |0,i,l) 



oil)' 

33] . The parameter Vq is fixed by requiring 



a G, 



30 
03 



cr 



0,0, |) 



The variance of this distribution is given by 



a 



MaxEnt 



Ggg(a 3 |- 1,0,0) G^( a 3 |0,i,l) 
KM 0,0,1)] 2 
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- 1. 



(12) 



(13) 



which tends to zero as a tends to oo (nearly chaotic regime) and behaves as a" 3 / 2 at small 
a (nearly integrable). Figure 1 compares the distributions of LDOS in Eq. (9) and (10) 
for same variances. The figure clearly suggests the equivalence of the two distributions are 
nearly equivalent, except at small values of a or lo that correspond to variances a 2 > (V) 2 , 
which is adapted in the superstatistical approach with nearly integrable systems. This should 
not make an obstacle for using the Altshuler-Prigodin LDOS distribution in superstatistical 
RMT. The superstatistical approach is not expected to model a system in the final stage of 
transition out of chaos, where the system approaches the integrability limit. As previously 
mentioned, this is a base-invariant approach. Integrable systems by definition have well 
defined complete set of eigenvalues that yields a diagonal Hamiltonian matrix when used as 
a basis. 



Schomerus et al [30| calculate the probability distribution of the local density of states 
in a disordered one dimensional conductor or single-mode waveguide. They start from the 
relation between the LDOS and the retarded Green function, which they have obtained by 
a numerical simulation for a wire with one or both fixed ends . They have given exact 
results for the distributions of the local densities of states in one-dimensional localization, 
contrasting the microscopic length scale (below the wavelength) and mesoscopic length scale 
(between the wavelength and the mean free path). Their data points are shown in Fig. 2 
for the case of a closed wire. These data result from a numerical simulation for a wire of 
length equal to 55 times the mean free path, with the LDOS computed halfway in the wire. 
Figure 2 also shows the result of calculation using the Altshuler-Prigodin distribution (9) by 
a dashed line and the distribution in Eq. (9) by a solid line. The mean LDOS is fixed to be 
equal 1 in the three distributions. While the Altshuler-Prigodin distribution fits the results 
of numerical simulation better, the superstatistics distribution still provides a reasonable 
representation of the data. 

The absence of quantitative agreement between the predictions of Eqs. (9) and (10) does 
not imply that the LDOS v is not equal to the inverse of local mean level spacing D. The 
Altshuler-Prigodin distribution (9) is derived from an elaborate theory that starts from the 
relation between the level density and the retarded Green's function for noninteracting elec- 
trons in a wire and arrives at a relation between the microscopic LDOS and the reflection 
coefficients, and then uses elaborate methods perform the local spatial average that gives 
the LDOS. It is interesting to note that the distribution in Eq. (10) is derived in lJ from 






S 



the principle of maximum entropy with the constraints of fixed (D~ 2 ) and (D) as mentioned 
above. If we replace the first constraint by another one that requires that (D^ 1 ), we ar- 
rive exactly to the Altshuler-Prigodin distribution. The fact that the Altshuler-Prigodin 
distribution fits the numerical-experimental data perfectly, as shown in Fig. 1, while the 
distribution (10) fits the data only qualitatively suggests that f(D) may better be obtained 
by maximizing entropy under the conditions of fixed (D) and (D^ 1 ). 



IV. NEAREST-NEIGHBOR-SPACING DISTRIBUTION 



The question that this section tries to answer is whether the distribution of LDOS ob- 
tained in the study of disordered metals is suitable for describing the transition out of 
chaos within the superstatistical approach to RMT. For this purpose we consider a system, 
which is described by a superposition of random matrix ensembles of different LDOS with 
probability-density function given by the Altshuler-Prigodin distribution Pap(uj, u ap ; v). Its 
spectral characteristics are expressed as averages of the corresponding characteristics of a 
Gaussian random ensemble in analogy with Eq. (5). Then, the superstatistical NNS distri- 
bution is given by 

/■oo 

p AP (u,s) = P A p{uj,u ap ;u) p WD (u,s) du, (14) 
Jo 

where Pwd(v, s) is a Wigner-Dyson distribution with mean level spacing equal to 1/u, 



Pwd(*A s) = ^z/ 2 sexp \-jP 2 s 2 ) . (15) 



We then obtain 



p(u,s) = -\[2nujv\ v s J \fx exp 
4 Jo 



CO ( 1 
00 — — I X H | r^AP^" 5 " 



dx (16) 



2 V xj 4 

the parameter z^ap is fixed bey requiring that the mean-level spacing is equal to 1 

which yields 

z/ap = 1 + 1/co. (17) 

Unfortunately, we are not able to evaluate the integral in Eq. (16) analytically, but there is 
no problem in its numerical evaluation. 

We shall demonstrate the quality of the NNS distribution in Eq. (16) by applying it 
to the results of a numerical experiment by Gu et al. [si| on a random binary network. 
Impurity bonds are employed to replace the bonds in an otherwise homogeneous network. 



The authors of Ref. 



31| numerically calculated more than 700 resonances for each sample. 



For each impurity concentration p, they considered 1000 samples with totally more than 700 
000 levels computed. Figure 3 shows their results for four values of concentration p. The 
figure also shows the best fits obtained for NNS superstatistical distributions obtained using 
the parameter distributions following the Altshuler-Prigodin formula and the ones obtained 
in [13J using the principle of MaxEnt. The high statistical significance of the data allows us 
to assume the advantage of the superstatistical distribution base on the Altshuler-Prigodin 
LDOS distribution for describing the results of this experiment. 



V. SUMMARY AND CONCLUSION 



Many characteristics of disordered metals and those of insulators have been understood 
using RMT. The spectral fluctuations of a disordered metal are well described by random 
matrix theory, while the fluctuations of a system in the insulating regime follow the Poisson 
statistics. The nature and the details of the metal-insulator transition, on the other hand, 
still belong to the most intensively studied problems. Localization of wave functions by dis- 
order can be seen in the fluctuations of the density of states, which can be probed using the 
tunnel resistance of a point contact. We have previously described the analogous transition 
between regular and chaotic dynamics within the framework of RMT as a superposition of 
two statistics, namely one described by the matrix-element distribution exp [— r/Tr [H'Hj^ 
and another one by the probability distribution f{jf) of the inverse variance of the matrix 
elements, which is derived from the principle of maximum entropy. In the present investiga- 
tion, we show an interesting connection between the distribution f{rj) and the distribution 
of LDOS which is used in the study of disordered systems. We have also found that the 
LDOS distribution that follows from f(j]) is in a reasonable agreement with the numerical 
numerical simulation of LDOS is a weakly-disordered metallic wire. The agreement is of 
course not as good as for the distributions obtained by Altshuler and Prigodin for LDOS in 
one-dimensional disordered systems using a more sophisticated approach. However, we ob- 
tain the Altshuler-Prigodin distribution if we modify the (optional) constraints imposed on 
entropy maximization in our previous work, by requiring a fixed value of {f] 1 ^ 2 ) rather than 
fixing (rj). Therefore, both the Altshuler-Prigodin distribution and the parameter distribu- 
tion used in our previous work are equally well based on the principle of maximum entropy, 
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both satisfy the Beck-Cohen criteria, but the comparison with experimental data prefers 
the former distribution. In conclusion, we suggest to use the Altshuler-Prigodin formula as 
a parameter distribution that describes the superposition of GOE's in the superstatistical 
ensemble. This has been illustrated by an analysis of a high-quality numerical experiment 
on the NNS distributions of resonance spectra of disordered binary networks. 
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Figure captions 

FIG. 1. (Color on line) Distributions of LDOS estimated by using the principle of MaxEnt 
and previously used in the superstatistical approach to RMT [13J (solid line) compared with 
the ones obtained by Altshuler and Prigodin for a closed one-dimensional wire 29] (dashed 
line) that have same variances. 

FIG. 2. ((Color on line) Distribution of LDOS in a disordered metal numerically calcu- 
lated by Schomerus et al. (histogram) compared with the corresponding distribution 
deduced from the superstatistical approach to RMT [13j (solid line) and that obtained by 
Altshuler and Prigodin for a closed one- dimensional wire [29] (dashed line). 

FIG. 3. (Color on line) Nearest neighbor spacing distributions of geometrical resonances 



in random network, calculated by Gu et al. 



3l| compared with the superstatistical distribu- 



tions in which the parameter distributions are estimated using the MaxEnt principle (solid 
lines) as well as those based on the Altshuler-Prigodin distribution of LDOS (dashed lines). 
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